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Tjj- ■ Abstract 

^ ■ The usefulness of the 5° ^ D-al,D- (pn' K+K-Ti',al p°7r+ ^ 

7r"'"7r tt"*" decay chain is investigated for the i^s-reconstruction in the future ATLAS 
■ iJ^-mixing experiment. It is shown that this decay channel is almost as suitable for 

\ this purpose as previously studied D^tt^. 

^ ; 1 Introduction 

, , . In the Standard Model 5° and 5° are not mass eigenstates. Instead we have (the small 

! CP-violating effects are neglected) 

„n -Bi + Bo -„ B^ — Bo , , 

B0 = ^ 8" = ^. (1) 

So the time evolution of the Bi states looks like 

Bi{t)=B,{{])exp^-'-{mi-i^)t^ , (2) 

where is the mass eigenvalue and Fj - the corresponding width. 

It follows from (1) and (2) that the probability for B^ meson not to change its flavour 
after a time t from the creation is 

P«»-°(i) = |(i3°(i)|B°(0)>P = ie-*(co*^ + cos^) . (3) 
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and the probability to convert into the 5° meson - 

P«'*(i) = l{S°W|B°(0)>P = ie-¥(co=h^-cos^) , (4) 

where F = ^(Fi + F2) is the average width and AF = Fi — F2 . So Am = mi — 1712 mass 
difference between the B mass eigenstates defines the oscillation frequency. Standard 
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Model predicts that ~ > 1, V^ij being the Cabibbo-Kobayashi-Maskawa 

matrix element. Therefore the mixing in the meson system proceeds much more faster 
than in the system. 

The total probability x that a B^ will oscillate into B^ is 

,./>n<(l-.^)4(^.^)(l-.^) .^^^..^f . (5) 

In the first S^-mixing experiments just this time integrated mixing probability was 
measured. The result = 0.69 ±0.07 shows that in the Bs system ^ 1 is expected. 
In fact the allowed range of Xg is estimated to be between ~ 12 and ~ 30 in the Standard 
Model 1^. Such a big value of Xg makes impossible time integrated measurements in the 
Bs system, because x i^i (5) saturates at ~ 0.5 for large values of x. 

Although it was thought that unlike the kaon system for the B mesons the decay width 
difference can be neglected 0, nowadays people is more inclined to believe the theoretical 
prediction that the b — > ccs transition, with final states common to both Bg and Bg, 
can generate about 20% difference in lifetimes of the short lived and long lived fi^-mesons 

i- 

But we can see from the (3 ~ 5) formulas that the effect of nonzero y is always ~ y"^ 
and so of the order of several percents, because ?/ ~ 0.1 is expected. In the following we 
will neglect this effect and will take y = 0, though in some formulas y is kept for reference 
reason. 

The development of high precision vertex detectors made it possible to measure ||^ in 
the Bd system the time dependent asymmetry 

dB^B" dB'^B'^ a j. 

COS — - — . (6) 



pBOBO _j_ pB°BO 

The same techniques can also be applied to the B^ — Bg system. 

Recently the ATLAS detector sensitivity to the Xg parameter was studied using 
5° — > D~Ti^ — > (f)TT~TT'^ —>■ K^K~7T~Tr^ dccay chain for Bg meson reconstruction. It 



was shown that Xg up to 40 should be within a reach |T0|. The signal statistics could be 
increased by using other decay channels, like B^ D^cit , B^g ~^ J /'^K*^. 

The purpose of this note is to study the usefulness of the decay chain B^ — > D ^ (i\ — > 
(j)7[~pTT'^ K^K^ 71^71^71^71^ for Bg mcsou reconstruction in the ATLAS i^s-mixing 
experiments. 

2 Event simulation 

About 20 000 following b-decays were generated using the PYTHIA Monte Carlo program 

Q 
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(p^ > 6 GeV/c, \r\< 2.2) fitag ^ hi ^ 5° ^ a+ 



The impact parameter was smeared using the following parameterized description of 
the impact parameter resolution 

aip = 14 © 72/(pT\/|sin^|) (x^ = 20 © 83/(pt\/| sin ^P) , (7) 

where resolutions are in /im and 9 is the angle with respect to the beam line. It was shown 
in that this parameterized resolution reasonably reproduces the results obtained by 
using the full simulation and reconstruction programs. 

For the transverse momentum resolution an usual expression 

= 5 . lo-V © 1.2% (8) 

Pt 

was assumed. 

Track reconstruction efficiencies for various particles were taken from []T0[. Because 



now we have 6 particles in the final state instead of 4 for the 5° — D~t[^ decay channel, 
we expect some loss in statistics due to track reconstruction inefficiencies, but the effect is 
not significant because the investigation in [|1^ indicates a high reconstruction efficiency 
of 0.95. 



3 Event reconstruction 

The topology of a considered 5° decay chain is shown schematically in a figure: 
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f^tag 



The Bs decay vertex reconstruction was done in the following three steps. 

First of all the D~ was reconstructed by finding three charged particles presumably 
originated from the D~ decay and fitting their tracks. For this goal all combinations 
of the properly charged particles were examined in the generated events, assuming that 
two of them are kaons and one is pion. The resulting invariant mass distribution is 
shown in Fig. la for signal events. The expected Dj peak is clearly seen along with 
moderate enough combinatorial background. Cuts on A(j)KK , ^^kk and \Mkk — 
were selected in order to optimize signal to background ratio. To select one more cut on 
\Mkkit — the information about the invariant mass resolution is desirable. Fig. 

2a shows the reconstructed D~ meson from its true decay products. The finite invariant 
mass resolution is due to applied track smearing and equals approximately to lOMeV/c^. 

After D~ meson reconstruction, af meson was searched in three particle combinations 
from the remaining charged particles, each particle in the combination being assumed 
to be a pion. Fig. lb shows a resulting invariant mass distribution for signal events. 
Because of huge width of af, signal to background separation is not so obvious in this 
case. If af is reconstructed from its true decay products as in Fig. 2b, its width is 
correctly reproduced. To draw out af from the background, further cuts were applied on 
Acf)^^ , Ae^^ , |M^^ - Mp\ and |M^^^ - Ma^\. 

At last decay vertex was fitted, using reconstructed and af. 

Almost the same resolution in the i?<j-decay proper time was reached ar ~ 0.064ps, as 
in 1^. The corresponding resolution in the B-meson decay length in the transverse plane 
is ~ 87/im. The relevant distributions are shown in Fig.3. 
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4 Signal and background 



Branching ratios and signal statistics for the 
Table 1. Note that we use an updated value for Bt{Dj - 



channel are summarized in 
branching 



■ 071") from [|T2 

ratios are still unknown experimentally. Neglecting SU(3) unitary symmetry breaking 
effects, we have taken Br(_B° D'af)^ Br (5° D^al). 

Table 1. 

Branching ratios and signal statistics for 5° D^al{12Q{)). 



Parameter 


Value 


Comment 


L [cm^'^s^^] 


1033 




t[s] 


10^ 




(j{hh)/a{tot) 


~ 1/100 




a{hl) 


~ 500 




a{hh /iX) 


~ 2.24 


pi^>Q GeV/c 






\r]^'\ < 2.2 


N{bb fiX) 


2.24 X 10^° 




Br{b 5°) 


0.1 




i?r(i?o ^ D;at) 


0.006 




Br{D- 4>7T-) 


0.035 




Br{(j) ^ K+K-) 


0.491 




Br{at p°7r+) 


~ 0.5 




5r(p° TT-TT+) 


~ 1 




N{K+K-TT-7l + 7l-7l + ) 


116000 





Acceptance and analysis cuts are summarized in Table 2. We take a track reconstruc- 
tion efficiency of 95% and a lepton identification efficiency of 80%, as in [10 . 
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Table 2. 

Analysis cuts and acceptance for D~ af (1260) (for 10^p6~^ integrated luminosity). 



Parameter 


Va ] u e 


domment 


AT( l<r+ T<"—'Tr — 'Tr + 'Tr — m-+\ 

iV i\ TT TT TT TT J 


i iOUUU 




Uuts . 






Pt > 1 Ge\//c 






|r/| < 2.5 






1\ \ / 




u.u /o 














1 A/f,^,, _ A/f,| ^9(1 A/fpT//r2 






\ -^^^^KKtv — -^^Z)~l IVlt-VfC 












A^^^ < 15° 






l/Vf _ /Vf nl < iqo MpVIc^ f+Sn-l 
TVf _ M + 1 < 800 MeVlc^ 






AT f 7^4- Ty— 

iV \K ^ K 71 TT^TT TT^J 


57o5 


5.0/0 


Vci Lc-X. ill ^ \ IZr.U 






a^^ vertex fit < 12.0 






5° vertex fit < 0.35 






proper decay time > 0.4 ps 






B^ impact parameter < 55 /im 
5° pt > 10.0 GeV/c 






N{K^K^ 71^71^71^71^) after cuts 


3505 


3.0% 


Lepton identification 


0.8 




Track efficiency 


(0.95)6 




N{K^K^7i^7i^7i^7i^) reconstructed 


2065 


1.8% 



As we see, about 2065 reconstructed 5° are expected after one year run at £ = 
^ luminosity. The corresponding number of events within one standard devia- 
22 MeV/c^) from the B^ mass equals 1407. This last number should be compared 



lO^^cm 



tion 



to 2650 signal events, as reported in |T^, when 5° D~7i^ decay channel is used. 

Events which pass the first level muon trigger [px > 6 GeV/c, \t]\ < 2.2) are predomi- 
nantly bb events. Background can come from other B decays of the same or higher charged 
multiplicity, and from random combinations with some (or all) particles originating not 
from a B decay (combinatorial background). 

The following channels were considered and no significant contributions were found to 
the background: 

• i?^ ^ D^af . These events don't pass the analysis cuts, because the D~ mass is 

shifted from the _D~ mass by about 100 MeV, and so does the B'^ mass compared 
?o 



to the B^ mass. 



Ab - 
from 



A; 
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'7T' followed by A+ ^ p/^^tt+tt+tt" . Taking Br{Ab A+vr") ^ 0.01 
, we see that the expected number of pKATi events, originated from this 
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source, is only five times less than the expected number of truly signal events. But 
the decay topology for this decay chain is drastically different (1+5, not 3+3) and 
therefore it is unexpected that significant amount of the B-decays will be simulated 
in this way. 

Note that even for Bg D^n^ decay channel the similar background is negli- 
gible 0, although Br{A^ pK^7c~^) is about 44 times bigger than i?r(A^ 

— s> D^af. About 10 000 such events were generated by PYTHIA and then 
analyzed. Using Br(i?2 — > D +ar) < 2.7 ■ 10-^ from p| and assuming that B^ — >■ 



tti decay goes through B^B^ oscillations: B^ ^ B^ ^ a^, and therefore 
Br{B^ D-at) = XdBr{B°^ Dfa{) < 4.3 ■ 10-^ we have got Fig.4. It is seen 
from this figure that because of Mb^ — Mb^ ~ lOOMeV^ mass shift, the contribution 
of this channel to the background proves to be negligible. 

Note that Fig.4 refers to the total number of the B^ D~af events. In fact the 
distribution of these events with regard to the decay proper time is oscillatory, Xd 
(not Xg) defining the oscillation frequency. So in general this will result in oscillatory 
dilution factor. The conclusion that this dilution factor is irrelevent relies on the 
fact that no candidate event was found with invariant mass within one standard 
deviation from the Bg mass for 6 ■ pb~^ integrated luminosity. 

A huge Monte-Carlo statistics is needed for combinatorial background studies. No 
candidate event with Mgo - IbOMeV/c^ < Mkktttttttt < Mbo + IbOMeV/c^ was found 
within ~ 3 ■ 10^ inclusive fiX events. This indicates that signal/background ratio is 
expected to be not worse than 1:1. 



5 Dilution factors 

The observation of the B — B oscillations is complicated by some dilution factors. First of 
all the decay proper time is measured with some accuracy a. From previous discussions 
we know that in our case a = 0.064ps is expected. Due to this finite time resolution, the 
observed oscillations are convolutions of the expressions (3) and (4) given above with a 
Gaussian distribution. For example 



pB^BO 




oo 

h 




^ 2 








oo 


1 _« 


cosh 


AF 








'2h 



An 



Ams 



cos ■ 



h 



exp 



it 



2(t2 



cr Am 

'-) - D time cos ^—[t 



a— 

t' 



ds 



2txo 



(9) 



where D 



time 



exp 



-) 



Hx" 



T = 7^ 



h 
r- 



So the main effect of this smearing is the reduction of the oscillation amplitude by 
Dtime- This is quite important in the Bg system where x ^ 1. There is also a time shift 
t — s> t — cr^ in (9). This time shift does not really effect the observability of the oscillations 
and we will neglect it. 
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In fact (9) is valid only for not too short decay times t ^ a, because in (3) and (4) 
distributions t > is assumed. 

Another reduction in the oscillation amplitude is caused by the particle/ antiparticle 
mistagging at t=0. In our case particle/antiparticle nature of the B meson is tagged by 
the lepton charge in the semileptonic decay of the associated beauty hadron. Mistagging 
is mainly due to 

• B — B oscillations: accompanying b-quark can be hadronized as a neutral B meson 
and oscillate into B before semileptonic decay. 

• b c cascade process, then the lepton is misidentified as having come directly 
from the i?-meson and associated to the 6 — >■ decay. 

• leptons coming from other decaying particles (K,7r,...). 

• detector error in the lepton charge identification. 

Let 7] be the mistagging probability. If we have tagged N B^ mesons, among them 
only {1 — ri)N are indeed B^-s and r]N are B^-s misidentified as B'^-s. So at the proper 
time t we would observe {P^°^\t) = P^°^°{t) due to CPT mvariance 



N 



;i-,)P-°-""(t)+,P-°-°(t)]=fe-x 



An Amt 

cosh — (1 — Zf]) cos — - — 

2h h 



decays associated to the B^ meson and therefore 



^ ' 2 



, Art ^ ^ Amt 
cosh (1 — 2r]) cos — - — 



So the dilution factor due to mistagging is Dtag = 1 — 21]. In our studies we have taken 
Dtag = 0.56, as in [T^. 

Finally the dilution can emerge from background. Suppose that apart from 

Nsignal -El f , AH Amt 

— e R cosh — cos — - — 

2 \ 2n h 

events with B ^ B oscillations we also have Nf,ack{t) additional background events. Half 
of them will simulate B meson and half of them B meson (assuming asymmetry free 
background). So the observed number of would he B ^ B oscillations will be 

Nsignal -Itf . ATt Amt\ Nhack{t) _rt( AVt Amt 

— e R cosh — cos — - — H ~ e r cosh — iJback cos — - — 

2 \ 21i h J 2 \ 2h h 

and the oscillation amplitude will be reduced by an amount 

Nsignal ■ cosh ^ 



D 



back 



Nsignal + Nbackit)^'' 



Neglecting the proper time dependence of this dilution factor (that is supposing that 
the background is mainly due to S-hadron decays and therefore has approximately the 
same proper time exponential decay as the signal [|l5]),we have taken D^ack ~ 0.71 which 
corresponds to the 2:1 signal/background ratio. 
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6 Prospects for Xs measurements 



For 6 ■ 10^ pb ^ integrated luminosity the number of reconstructed -B^-s would reach 
~ 8000 from the analyzed channel alone. Another ~ 16000 B^-s are expected from the 
B° D-TT+ channel |, |ig. 

For events in which meson does not oscillate before its decay, the Dg meson and 
the tagging muon have equal sign charges. If the -8° meson oscillates, opposite charge 
combination emerges. The corresponding decay time distributions are 

dn{++) N f Xst. 
— — e ^1 + 19 cos — 



dt 2r V r 

v^ri-Dcos(i^)) (10) 



dt 2t 

D is the product of all dilution factors and is the total number of reconstructed -B^-s. 

The unification of samples from B^ D~af and B^ decay channels allows 

to increase Xs measurement precision. 

Fig. 7 and Fig. 8 show the corresponding 



dn{++) dn{+-) 



, X at , 



dt ~'~ dt 



asymmetry plots for Xs = 20 and 35. 



7 Conclusions 

It seems to us that B^ ^^(^t decay channel is almost as good for the i?s-mixing 
exploration as previously studied 5° D^tx^ and enables us to increase signal statistics 
about 1.5 times. Further gain in signal statistics can be reached |10| by using 5° 
J/i/jK* decay mode and considering other decay channels of D^. These possibilities are 
under study. 

We refrain from giving any particular value of attainable upper limit. Too 

many uncertainties are left before a real experiment will start. Note, for example, that 
about two times bigger branching ratios for both Bg D^n^ and Bs D^al decay 
channels are predicted in [|T6|. ~ 500/i6 as a hh production cross section can also have 
significant variation in real life ||l^ . 



So although the results of this investigation strengthen confidence in reaching Xs as 
high as 40 [|1^, it should be realized that some theoretical predictions about iJ^-physics 



and collider operation were involved and according to T.D.Lee's first law of physicist JTS] 



"without experimentalist, theorist tend to drift". However maybe it is worthwhile to 
recall his second law also "without theorist, experimentalists tend to falter". 
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Figure 1: Invariant mass distributions of three charged particle combinations in signal 
events, assuming 2K + tt (a) or Stt combination (b) as described in the text. 
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Figure 2: Three particle invariant mass distributions of reconstructed and events. 
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Figure 3: Proper time (a) and transverse radius (b) resolutions for the reconstructed B[ 
decay vertex. 
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Figure 4: Six particle invariant mass distribution corresponding to the meson. Dashed 
hne - expected upper hmit for background from S° decay. 
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Figure 5: Asymmetry distributions for B'^ — > (a), — * 7r+ (b) and when both 

channels are used (c), for Xg — 20. 
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Figure 6: Asymmetry distributions for B'^ — > (a), — * 7r+ (b) and when both 

channels are used (c), for Xs — S5. 
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